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ABSTRACT 

Even highly improved variants of lattice QCD with staggered fermions show 
significant violations of taste symmetry at currently accessible lattice spacings. 
In addition, the "rooting trick" is used in order to simulate with the correct 
number of light sea quarks, and this makes the lattice theory nonlocal, even 
though there is good reason to believe that the continuum limit is in the correct 
universality class. In order to understand scaling violations, it is thus necessary 
to extend the construction of the Symanzik effective theory to include rooted 
staggered fermions. We show how this can be done, starting from a generalization 
of the renormalization-group approach to rooted staggered fermions recently 
developed by one of us. We then explain how the chiral effective theory follows 
from the Symanzik action, and show that it leads to "rooted" staggered chiral 
perturbation theory as the correct chiral theory for QCD with rooted staggered 
fermions. We thus establish a direct link between the renormalization-group 
based arguments for the correctness of the continuum limit and the success of 
rooted staggered chiral perturbation theory in fitting numerical results obtained 
with the rooting trick. In order to develop our argument, we need to assume the 
existence of a standard partially-quenched chiral effective theory for any local 
partially-quenched theory. Other technical, but standard, assumptions are also 
required. 



I. INTRODUCTION 



On a hypercubic lattice in four dimensions, the continuum limit of lattice QCD with 
staggered fermions [l| contains four "tastes" of mass-degenerate quarks per staggered fermion 
field 0,1,1, 9-^ Hence, if we introduce a separate staggered fermion field for each physical 
light-quark flavor (up, down, and strange), the continuum limit consists of QCD containing 
four up, four down, and four strange quarks. 

A simple solution to this problem is to adjust for the excessive multiplicity by taking the 
fourth root of the fermion determinant for each staggered fermion field Heuristically, if 
the staggered determinant factorizes into four identical determinants in the continuum limit, 
one for each taste, taking the fourth root corrects for the taste multiplicity. The desired 
theory, QCD with one up, down and strange quark each is then obtained in the continuum 
limit. Since the staggered determinant is positive for any real, nonzero bare quark mass m, 
and the continuum determinant is (formally) positive for positive quark mass, the positive 
fourth root should be chosen.^ The continuum quark mass is proportional to |m|, which 
undergoes only a multiplicative renormalization, because staggered fermions have one exact 
chiral symmetry. 

This procedure, the "fourth-root trick," raises a number of questions jo], [l^, 11]. The 



fourth root of a determinant cannot in general be written as a Grassmann integral with a 
local action. Therefore, the first question is whether the theory defined by the fourth-root 
trick is local and unitary. 



In Ref. [I4] we showed that, as might be expected, the fourth-root staggered theory is not 



local at nonzero lattice spacing a. Continuing correlation functions defined in the Euclidean 
theory to Minkowski space will lead to violations of unitarity at a 7^ 0, on a distance scale 
set b y t he lightest particles in the theory, the Goldstone bosons. For examples of this, see 
Ref. [ill, as well as Sec. 6 of Ref. jl4|, which we will revisit later in this paper. 

The origin of these diseases can be traced back to the taste symmetry-breaking part of the 
staggered Dirac operator. This taste-breaking part corresponds to a dimension-five irrelevant 
operator. Thus, in the local, unrooted staggered theory, all taste symmetry-breaking effects 
are expected to vanish in the continuum limit, where exact U{4) taste symmetry will be 
restored for each of the four up, four down, and four strange quarks present in that theory. 

The leading power-law scaling of irrelevant operators is characteristic of any local and 
renormalizable theory, such as in particular the unrooted staggered theory. This brings 
us to the second question: Does the same scaling persist in the fourth-root theory? Two 
related considerations make it natural to address this question via a Renormalization-Group 
(RG) approach. To begin with, the RG framework allows us to define what we mean by the 
continuum limit. This is done by performing n + 1 blocking steps^ on the original lattice 
theory, with its fine spacing a/ = a, each time increasing the lattice spacing by a factor 
of two, to arrive at an RG-blocked theory formulated on a lattice with a coarse spacing 
ttc = 2"+^aj. Keeping Uc fixed and small in physical units, ^ ^qcd^ while sending 
n —>■ 00 (and thus af 0), one obtains a coarse-lattice theory describing the continuum 
physics. An RG framework is also natural because the restoration of taste symmetry is 
only expected to occur on distance scales much larger than the original lattice cutoff aj. RG 



^ We assume the usual choice of only a single-site bare mass term. 

^ For the case of an odd number of quarks with negative quark rnass, see Refs. [Tj, |8|. 

^ See Sec. Illll for an explanation of the convention adaf = 2"+-'^ [15|. 
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blocking removes the short- distance fluctuations while modifying the action of the remaining 
degrees of freedom by local terms only. When we increase the number of blocking steps n, 
the blocked theory becomes more taste symmetric, and we eventually recover exact taste 
symmetry in the continuum limit n — oo. 



Using this RG framework, it was argued in Ref. [15|| that the continuum limit of QCD 
with rooted staggered fermions is a local theory that belongs to the correct universality class. 
There are strong arguments that the fourth-root theory, while nonlocal, is nevertheless 



renormalizable |10l . Il6l . Il7l |. and this is the fundamental reason behind the validity of its 
continuum limit. The detailed reasoning is based on a number of technical assumptions, all 
of which are very similar to the assumptions needed to establish the nature of the continuum 
limit for the unrooted staggered theory. Further analytic and numerical work aimed at 
confirming the technical assumptions of Ref. p^] would add direct and strong evidence for 



the validity of the fourth-root trick. For full details, we refer to Ref. [15 



intuitive accounts, we refer to Refs. 
scaling of the taste-breaking effects 



m 21 



^ ^ for shorter, more 

lol . 11|. We stress that one key element — the anticipated 
-has been corroborated by extensive numerical studies 



Assuming that the rooted staggered theory has the correct continuum limit, this leaves 
us with a third question. While the anticipated scaling of taste-breaking effects is observed, 
these effects are clearly not negligible at present [l3, 18, 19, 2^, 21|. It is therefore imperative 
to take lattice artifacts into account in the effective continuum field theories (EFTs) such 
as the Symanzik effective theory (SET) or chiral perturbation theory (ChPT). The latter 
provides a central tool for analyzing the numerical data and performing the chiral and 
continuum extrapolations in the light-quark sector. In the case of rooted staggered fermions, 
we thus need to construct EFTs that take the discretization effects into account, including 
those that correspond to the nonlocal behavior of the theory at a 7^ 0. The construction of 
such EFTs is the subject of this paper. 

For the pseudo-scalar Goldstone-boson physics, a candidate EFT already exists; it is 
provided by staggered ChPT 2^ with the replica rule (rSChPT), or "rooted staggered 
ChPT" [23| . (Extensions to higher order and to heavy-light meson j25| and baryon [2^ 
rSChPT were recently given.) An argument for the validity of rSChPT was presented in 



Ref. [13], and reviewed in Refs. 0, El. The key feature of Ref. [14 1 is that the argument 
takes place completely within the context of chiral effective theories, and the replica rule is 
justified only in that context. Here we will need to introduce a somewhat different version 
of the replica rule, which will be justified in addition at the level of the fundamental lattice 
theory, but which will ultimately give the same results in the chiral theory. A detailed 
comparison of the two approaches will be made in Sec. IV C[ 

The overall goal in the current paper is to extend the standard procedure for the con- 
struction of ChPT for a local lattice theory to QCD with rooted staggered fermions. The 
standard procedure consists of two steps. The SET [l^] is constructed first. This can be 
done order by order in perturbation theory, but it is generally assumed that the SET is valid 
nonperturbatively as well. We will assume throughout that this includes partially quenched 
theories [l6|. In particular, we will assume that locality suffices, and that unitarity (which 
may be lost in partially quenched theories) is not necessary. Once the correct form of the 
SET has been established, its symmetries can be used to construct ChPT. Since the SET 
organizes the low-energy effective theory as a systematic expansion in the lattice spacing, 
one automatically obtains the chiral theory as an expansion in the lattice spacing as well. 

Establishing that EFTs can be constructed following the usual rules for QCD with rooted 
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staggered fermions thus constitutes a fundamental step in understanding the effects of root- 
ing at nonvanishing lattice spacing. The main thrust of this paper is the construction of the 
SET for the rooted theory; obtaining the corresponding ChPT is then straightforward, and 
we show that it is indeed given by rSChPT. We emphasize that our construction applies 
to all commonly used versions of staggered fermions: standard (unimproved) staggered [H, 
Asqtad [l^l, HYP [2^, Fat7bar [sO?], HISQ [sT], etc. The only requirement is that the action 
have the usual staggered symmetries. The size of the discretization effects is of course dif- 
ferent with different versions of staggered fermions, but their form (and appearance at each 
order in a/) is the same. 

It is also important to note that the effective theories we ultimately construct are those 
for the relevant rooted staggered theory on the original (fine) lattice. The RG framework 
is used only as a tool in the derivation of these effective theories. Nevertheless, it is an 
indispensable tool: the conclusions of Ref. [l^ have to be valid in order for our construction 
of the EFTs to make any sense. We will assume this to be the case. 

The difficulty in constructing EFTs for the rooted theory is the following. Consider for 
simplicity a staggered theory with a common power, denoted n^, of the fermion determinant 
for each staggered flavor in the theory. As long as rir is a positive integer the lattice theory is 
local, and the construction of EFTs proceeds as usual. In order to arrive at the fourth-root 
theory,^ however, we must set = 1/4. Our task is to ensure that a replica continuation 
may be performed: a well-defined procedure must be devised to reach the value = 1/4 at 
the level of an EFT, and the procedure must be consistent with the n^-dependence of the 
underlying lattice theory. 

In a diagrammatic EFT calculation, the dependence on the number of (sea) quarks arises 
in two ways. First, there is explicit dependence arising through loop diagrams. In addition, 
the coupling constants of the EFT (the Symanzik coefficients in the case of the SET, and 
the low-energy constants in the case of ChPT) depend in an unknown way on the underlying 
lattice theory, including in particular on the number of replicas Uj.. It is the latter dependence 
that makes our task nontrivial. In principle, one may envisage two basic obstructions to 
the replica continuation of the coupling constants in the EFT. Mathematically, a unique 
analytic continuation off the positive integers (which in the case at hand is where the theory 
is local) does not exist. Also, it could be that the replica continuation we have in mind will 
encounter a singularity precisely at the desired point = 1/4. 

The dependence of the underlying lattice theory on the number of replicas is both 
perturbative and nonperturbative; this means that proving that no obstacle to the replica 
continuation is present would be tantamount to solving the theory nonperturbatively. The 
key observation that makes our task nevertheless tractable is that, after a large number 
n of RG blocking steps, the taste-symmetry breaking effects are very small: the unrooted 
staggered theory with integer is very close to a f/ (4) taste-invariant theory. The rooted 
theory, with = 1/4, is then also very close to a local lattice theory, for which the standard 
construction of EFTs is valid. Indeed, the "re-weighted" taste-invariant theories introduced 
in Ref. [l5| are local whenever is a multiple of 1/4. The proximity of these local theories 
makes it possible to construct the SET and, later, ChPT, for the rooted theory. 

We will reach the SET for the rooted theory starting from the SET for the corresponding 
re-weighted, taste-invariant theory. The flavors of the taste-invariant theory will always 



^ The discussion generalizes easily to the isospin limit = = nig, where one takes the square root of 
a single staggered flavor with (bare) mass mi. 
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be kept in one-to-one correspondence with those of the continuum-hmit theory. In the 
taste-invariant theory the dependence of the Symanzik coefficients on the physical quarks 
is nonperturbative, and unknown, as usual. This does not pose any difficulty, because the 
number of physical flavors is never varied. 

During the intermediate steps of the derivation the parameter will take on a related, 
but different technical meaning. The precise definitions will be given and explained in Sec. lIIII 
below. As already mentioned above, we first approximate the staggered theory by a local, 
taste-invariant theory that belongs to the correct universality class. The (rooted) staggered 
theory will then be reached from the taste-invariant theory by "turning on" smoothly the 
taste-breaking effects. The dependence on Ur of the lattice theory will come only from 
the taste-breaking effects, which are nonlocal (for noninteger rir) but small. The difference 
between corresponding taste-invariant and staggered theories is of order the fine lattice 
spacing a/ of the original (unblocked) lattice. This will allow us to show that all the lattice 
correlation functions are polynomials in to any fixed order in the expansion in a/. The 
degree of the n^-polynomial is less than the order of the a^-expansion. The ra^-dependence of 
the Symanzik coefficients can then be determined unambiguously. It follows that the replica 
continuation is nowhere singular in the complex Ur plane, to the given order in a/. Finally, 
after performing the replica continuation, the parameter resumes its original role as the 
power of the staggered determinant in the lattice theory. The further transition to ChPT 
is essentially a repeat of the same reasoning. As will become clear later on, we do have to 
assume that a chiral effective theory can be constructed for any local, but partially-quenched, 
theory. This was already emphasized in Refs. [l^, 14 . 

The outline of this paper is as follows. In Sec. [TTlwe consider the symmetries of staggered 
fermions in some detail. We derive the form in which shift symmetry is realized in 
the SET, and thus in any other EFT derived from the SET. A quick overview of the most 
important observations of that section is given at its beginning, and any reader not interested 
in the details can skip the remainder of the section. 

In Sec. mil we come to the main part of this paper, the construction of the SET for 
QCD with rooted staggered fermions. We generalize the staggered theory to a class of 
partially-quenched theories in which it is possible to implement the program outlined above. 
In Sec. HV] we discuss the SET to quadratic order in the lattice spacing in more detail, in 
order to illustrate the general construction. In Sec. |V] we make the transition to the chiral 
effective theory, and demonstrate that it is indeed given by rSChPT. As an example, we work 
out in rSChPT the leading-order contribution to the connected scalar two-point function. 



following the calculation in Ref. [ij]. We then compare the present derivation of rSChPT to 
that of Ref. [13], using the respective discussions of the scalar two-point function to make 
the comparison concrete. The final section contains our conclusions. A brief account of this 
work was presented at Lattice 2007 33 • 



II. SYMMETRIES OF THE SYMANZIK EFFECTIVE ACTION FOR STAG- 
GERED FERMIONS 

Here we discuss the symmetries of unrooted staggered fermions that are most relevant 
for this paper, and the way they appear at the level of the SET. We begin with an overview 
of the main results of this section. In the following subsections we will then give a more 
detailed discussion. 
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1. The staggered fermion action is invariant under shift symmetry, which, in the contin- 
uum hmit, enlarges to the product of SU{4:) taste symmetry and translation symmetry. 
At the level of the SET, the taste part of shift symmetry takes the form of the 32- 
element group r4 generated by a set of four-dimensional Dirac gamma matrices 
with 

{e^,U = 25^. , /i,i/e {1,2,3,4}. (2.1) 



This result was derived to order in Ref. 22|. Here we give a general argument that 



makes it clear that the result is true to all orders in a. On the continuum quark fields 
q used in the SET, the generating elements of r4 can be chosen to act according to 

q^^fifi , q-^liti ■ (2.2) 

Here the field qjsh has a Dirac spin index (3 and an S'?7(4) taste index 6, with the 
matrices acting on the latter. 

On the lattice, a taste-basis field if) carrying the same indices as the continuura quark 
field q is related to the one-component field x by a unitary transformation [s], [33| 

^ = Qx, V^ = xQ^ (2.3) 

The field lives on a coarse lattice whose spacing is twice that of the original stag- 
gered action. The ultra-local, unitary matrix Q maps the one-component variables x 
on the sixteen sites of each even hypercube to the sixteen components of ip on the 
single corresponding coarse-lattice site. The transformation Q is required to be gauge 
covariant, and its choice is not unique. As a result hypercubic rotational symmetry 
is somewhat complicated in the taste basis. ^ Of course, since the one-component and 
taste bases are related by a unitarity transformation, the physical consequences of all 
staggered symmetries are preserved. 

A somewhat different taste-basis operator, that we will refer to as the "RG taste-basis" 
Dirac operator |34|], is defined by a Gaussian smearing of the unitary transformation 



(12.31) . The resulting inverse Dirac operator satisfies 

AaL = - + QD:t\Q^ , (2.4) 

where D^tag is the Dirac operator in the one-component formulation, and a is a param- 
eter of order 1/a. Even though the theories described by D^tag and Dtaste are no longer 
related by a simple, unitary basis transformation, they are physically completely equiv- 
alent, because the propagators differ only by a contact term. The advantage of the 
Gaussian-smeared transformation is that discarding the taste-breaking part of Dtaste 



does not introduce any fermion doublers [12, l3J] . Because the staggered theory and 
the taste- invariant theory have a similar fermion content, one can interpolate smoothly 
between them. This will prove useful for the derivation of the SET. 

In the one-component formulation, shift symmetry is a unitary transformation on the 
fields X and x {cf. Eq. (12. lip in the next subsection). Since Q is unitary, the same is 
also true for the fields and tp, and from this it follows that the theory in the RG 
taste basis is also invariant under shift symmetry. 



^ For a detailed discussion of rotational symmetry in this framework, see Ref. [l5| . 
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3. Because of staggered symmetries, discretization errors for theories with staggered 
fermions start at order a} |35|. This is not obvious if one considers staggered fermions 



in the taste basis of Refs. 
of Refs. 



15, 34 



or in the modified form used in the RG analysis 
where taste-breaking terms occur in the action starting at order a. 
In this case, shift symmetry connects the leading, taste-invariant term in the lattice 
action with the order a taste-breaking term, z.e., their relative strength is fixed. There 
exists a local field redefinition that brings the taste-basis lattice action into a form 
where the taste violations are explicitly of order a^, and shift symmetry is realized as 
in Eq. (12.21) . again up to order terms [36|. More generally, the momentum-space 
basis used in the derivation of Eq. (12.21) can be related to the taste basis by a non-local 
field redefinition. Because the construction of the SET proceeds order by order in a, 
the field redefinition in effect becomes local. Therefore, the SETs constructed in the 
taste basis and in the staggered (or momentum-space) basis are always related by a 
local field redefinition. 



4. Staggered fermions have an exact chiral symmetry when m 
f/(l)e symmetry, taking the form 0] 



0, often referred to as 



eix 



\Xi+X2+X3+X4 



(2.5) 



For 171 = 0, this implies that 



{Dtaste, 75 ® ^5} = -Aaste (75 ® ^5) Aaste 

a 



(2.6) 



where 75 acts on the spin index, and ^5 = ^i^2^3^4_acts on the taste index [12|, |3J]. In 
other words, Dtaste is a Ginsparg- Wilson operator [37| with respect to f/(l)e symmetry. 



Before we proceed, we return to the relation of our analysis and that of Ref. 2^- The 
SET at order was determined in Ref. 22j by enumerating the allowed dimension-6 lattice 
operators consistent with the lattice symmetries, including shift symmetry. It was then 
shown that shift symmetry is represented on the corresponding continuum operators as a 
r4 symmetry. A more direct method of determining the SET, which we follow here, is to 
enumerate continuum operators. This leads to the result of point 1, that shift symmetry 
always implies a taste r4 symmetry of the SET. 

In the subsections following below, we will discuss some of these observations in more 
detail. These subsections are not needed for the construction of the SET for rooted staggered 
fermions, which can be found in Sec. IIIII 



A. Diagrammatic argument 

Our first argument for claim 1 above is essentially perturbative, and assumes that we 
are working in the momentum-space representation of the one-component basis. This result 
may be considered a corollary of Ref. . To keep it self-contained, a summary of relevant 
facts from Ref. [3] has been included in the discussion below. 

We will consider diagrams with n external fermion and r external gauge-field lines, corre- 
sponding to an operator which appears at a certain order in the SET. On the lattice, because 
of the phase factors which appear in the staggered action, momentum is conserved modulo 
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TT (in this section we work in lattice units), and any such diagram will have an overall delta 
function for momentum conservation of the form 

+ ■ ■ ■ + + A;i + • ■ ■ + A^r + n) , (2.7) 

where 11 is a vector with components or vr. The delta function is the periodic delta function 
with period 2it. The (lattice) quark and anti-quark momenta are pi, i = 1, . . . ,n and the 
gluon momenta kj , j = 1, . . . ,r. 

Because we are interested in an operator in the SET, we may take all physical external 
momenta small. Fermion doubling then implies that on every quark line we need to split 
the momenta as 

Pi = qi + 7iA, , (2.8) 
in which qi lives in the reduced Brillouin zone (— 7r/2 < g^^ < n/2), and vta^ = vrAj, with 

A, e {(0, 0, 0, 0), (1, 0, 0, 0), ... , (1, 1, 1,1)}. (2.9) 

We now take all physical momenta and kj small — so small that their sum has no 
components as large as ±7r. The delta function in Eq. (12.71) thus factorizes into 

S{qi + --- + qn + k, + --- + kr) 6{7CA, + --- + 7CA„+Il) . (2.10) 

Now consider what happens to this diagram under a shift 

X{x) C^{x)x{x + fl) , (2.11) 

^^{X) = (-l)-M+l + - + -4 _ gi-CM- , 

where the last equality defines vr^^. In momentum space (with x{x) = jpG^^^xip))^ this 
takes the form 

X{P^) = X{q^ + vtaJ ^ e'('''+'^^-)"x(g. + vr^. + H.) ■ (2-12) 

Applying a shift in the fi direction to all external legs of our diagram, and noting that 
the jth external gluon line is multiplied by a factor e*^^^^'' under a shift, we obtain the total 
factor 

gj(giH hgn+fciH hkr)^ ^2 X3^ 

which, by virtue of the first delta function in Eq. ( I2.10p . is equal to one. Therefore, we 
may omit these (small-momentum) phase factors in the shift (12.121) . We conclude that the 
diagram is invariant under the modified symmetry 

xiQi + ^aJ^ e'^^'^^^-xiqi + VTA, + TTcJ , 2 = 1, . . . , n , (2.14) 

which does not act on the gluon fields. The transformation fl2.14p generates a representation 
of the group r4 acting on the quark fields. Indeed, applying the transformation first in the 
fi direction, and then in the u direction, one obtains (dropping the index i) 

X{q + tta) ^ e^^^-^+'^c-)'^ e'^^^'^^xiq + tta + vrc, + vr^ J . (2.15) 



8 



For /i = z/, we have (vr^^)/^ = (c/. Eq. (12. lip ), and Eq. fl2.15p thus reduces to the identity. 
For /i 7^ i^, 



C^{x + u) = C^{x) c'^^^m)- = +1 , ^ > zy , (2.16) 
Ux + u) = -Ux) ^ e'^''^^ = -1 , /i < z/ , 

implying that shifts anti-commute, just hke the generators of We may make contact 
with Eq. (12.21) by introducing 

MQ) = x{q + 'rfA) . (2.17) 
The transformation Eq. (12.141) can now be written as 

Mq) ^ $^(s^)AB0B(g) , (2.18) 

B 

for some 16 x 16 matrices satisfying the Dirac algebra 

{S^,S4 = 25^, . (2.19) 

Finally, we can perform a basis transformation such that = 1 (S) and transform back 
to position space to obtain Eq. (12. 2p . 

Our argument shows that the diagram is invariant under the symmetry (12. 2p if it is 
invariant under shift symmetry (12.120 . The group r4 may thus be used to restrict the form 
of the SET in accordance with the shift symmetry of the underlying lattice theory. This 
is a considerable simplification, because the group r4 does not mix operators of different 
dimensions, i.e., of different orders in the Symanzik expansion. 

The same reasoning goes through in a theory in which the staggered fermion fields carry a 
flavor index i = 1, . . . , n/: one simply labels the fields Xi and Xi in Eq. (12.110 with the extra 
index i. Since the gauge fields also transform under shift symmetry, the same shift symmetry 
acts on all staggered fields simultaneously. It thus follows that the discrete symmetry r4 
acts in the same way on all staggered fields Xi^ and does not enlarge to the group (r4)"^ 

As an aside, we note that the invariance of the diagram under shift symmetry has im- 
plications for the second delta function in Eq. (12.100 . Naively, it would seem to follow that 
n just has to be equal to the sum over all ha,, but in general this is not sufficient. The 
reason is that the vertex can contain explicit periodic functions of the external momenta, 
which leads to additional sign factors under a shift. This is best illustrated with an example. 
Consider a lattice vertex of the form 

J2 '^(^l + ^2 + k)6{7TA + 7TB + 11) COS {Qi + k + 'KA)vX{q2 + B)x{qi + '^A)A^{k) , (2.20) 
A,B 

in which we split pi = gi + vr^, P2 = ?2 + tt^, and take gi,2 and k to be small. Performing a 
shift on the x and x fields results in (dropping a factor 5{qi + q2 + k)) 

5{tta + vtb + n) cos {qi + k + tta)u e'("^+"s)'' (2.21) 

A,B 

X x{q2 + + 7r^^)x(gi + t^a + vr^JA^(/i;) 

= 5{tia + ttb + n) cos {qi + k + txa)v e'^'^^^" e*"" x{q2 + T^B)x{q\ + T^A)A^{k) , 
a,b 
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where we used that (vtj^^J^ = and that 271^^ = mod 27r. The vertex is thus invariant if 
^fj, + (^Cm)i' ~ ^ laiod 271. An example of such a 11 is vr^^, which is defined by the phase 
factors which appear in the staggered action: 

r]^{x) = e^"""^ = . (2.22) 



B. Group-theoretical argument 

There is a very simple group-theoretical way to derive the same result. Let 5^ be the 
shift in the fi direction. All elements of the shift-symmetry group can be generated from 
the basic four shifts, and it is thus sufficient to consider only the 5*^. In any irreducible 
representation of the group, 5*^ looks like 

S, ^ e'^-E^ , (2.23) 

with — 7r/2 < < 7i/2 the physical momentum in lattice units, and the matrices generate 
a representation of r4 [38|. All irreducible representations are either "bosonic," if each is 
mapped onto ±1 (sixteen choices), or "fermionic," if the are chosen to satisfy the Dirac 
algebra (12.191) . Any field appearing in an EFT for the staggered theory (such as the SET or 
ChPT) transforms in some representation of S*^ under a shift {i.e., with some choice of 
and S^). 

Now we use that any continuum EFT is also invariant under continuum translations, 
which, on a continuum field $ with momentum q, act as 

$(g) ^ e^'^-'^iq) , (2.24) 

for a translation over a displacement r. We may thus choose r such that q ■ r = — g^, follow 
by this translation, and again obtain a symmetry of the EFT. This symmetry is precisely 
the one generated by the S^, i.e., a representation of 



Taste basis 



The arguments in the previous subsections made use of the momentum basis of the one- 
component formalism. The Feynman rules for the staggered theory in the one-component 
basis j3] were (assumed to have been) used in the derivation of the SET. Also, the group- 
theoretical argument works naturally on the momentum basis, since that is where irreducible 



representations of the staggered symmetry group live [38|. Alternatively, one could have 
started from the taste basis. The SET derived from the taste basis will not look the same 
as that derived from the one-component formalism; but the two SETs should be physically 
equivalent. Since the one-component and taste bases are related by a (nonlocal) unitary 
transformation in momentum space [5], one expects that the SETs derived from them, too, 
will be related by a field redefinition. Moreover, to any finite order in a, the SET-level field 
redefinition should be local, because the same is true for the unitary transformation between 
the two bases, when expanded to the corresponding finite order in a. 

We illustrate this in the free massiess theory, working to order in the Symanzik ex- 
pansion. On the taste basis, shift symmetry takes on the form js], [Sg] 

^{y) ^ I ((^M + l5l^.^5)Hy) + (^M - 757m^5)V^(i/ + fi)) . (2.25) 
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The field ip, introduced in Eq. (12.31) . is in this case given exphcitly by 

^Pb{y) = ^ 5^(7a);36X(2i/ + A) , (2.26) 

A 

where 7^ = 7i^^7^^7^^74^''; and A runs over the set (12.91) . The normahzation factor in 
Eq. (12.261) differs from that in Ref. |[5| because we take ip to be in lattice units of the coarser 
lattice; whereas Ref. jSj] works in physical units. In momentum space, Eq. (12.25^ looks like 

^(p) ^ e^p./2 (^^ cos - 2757^^5 sin {p^/2)) ^{p) (2.27) 



= e'"'" [i, - -i,i,i,p, + 0{p')j iPip) . 

The factor e*^'"'^ corresponds to the factor e*'''' in Eq. (12.121) . because the lattice spacings 
differ by a factor two. Dropping the factor e*^**/^ on the same grounds as in Sec. Ill Al it is 
easily verified that the transformation (I2.27P becomes a generating element of and that 
it is a symmetry of the order-a SET in the taste representation, 

Sfree = f ^^P) (^^^^^ + \ ^^^^^^^1 + ^^P'^) ^(P) " (2.28) 

This may also be written as 

^ [ii [V^(p)z7mPm^(p)] + 2 [V^(P)75P'V'(P)(757/.)1 + 0{p^)\ , 

(2.29) 

where we consider the field ip^h as a 4 x 4 matrix. 

In momentuni space, the transformation relating the one-component and taste represen- 
tations is [i, [3^, 33 

^(P) = ^E(-l)^'''^A0B(?)e- 



A,B 

-iq-A 

. 1A'VB\H)'- 

A,B 



where again A and B take values in the set (12. 9p . and where q = p/2. The field (/)(g) was 
defined in Eq. (I2.17p . The transformation (I2.30p is indeed nonlocal, but its expansion to 
any finite order in a is local. For instance, upon expanding e^*"^'^ = 1 ± iq ■ A + 0{q^), and 
starting from Eq. (I2.29p . this field redefinition brings the action (12.280 into the form 



Sfree = E E 7^ '^aI?) {"^{"^ ,) AB + 0{q')) Ml) , (2-31) 



H AB 



where the T ^ matrices form a 16-dimensional representation of the Dirac algebra and com- 
mute with the taste matrices Sj, defined in Eq. (12.180 . Note that Eq. (I2.3ip is expressed in 
units of the fine lattice spacing. 
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Let us also briefly consider the RG taste representation deflned by Eq. (12.41) in the massless 
free theory. To order the action is given by [15] 

E f ^(p) (w. 75^5^.^; + o{p')^ m . (2.32) 



This action is invariant under U{1)^ symmetry in "Ginsparg-Wilson-Liischer" (GWL) [37|, 



40[ | form. In the free theory, this symmetry looks like (again to order a) 



5^{P) = 75^5 I 1 - - 5Z + (^(P') 1 ^(P) ' (2-33) 



a 



5^/,(p) = 7/'(p) 75^5 . 

In this case, we may flrst carry out a fleld redeflnition 



1 + , (2-34) 

M / 



^(p) i){p) , 

followed by (12.301) . to bring the action into a form without terms of order a. Note that 
Eq. (I2.34P is nothing but the free-fleld, order-a form of the fleld redeflnition 

{I- D/aY^i) (2.35) 

with here D = Dtaste, which transforms the GWL form of U{l)e symmetry into a standard 
75^5 symmetry 0|. 

As a flnal note, we observe that to this order in a, fleld redeflnitions can be carried out 
such that the resulting action is invariant under the full U{4) taste symmetry. This turns 
out to be true to all orders in a in the free theory |^], but not in the interacting theory. 



III. DERIVATION OF THE SYMANZIK EFFECTIVE ACTION 

We begin by considering a theory with replicas of one staggered fermion with bare 
mass m, in the RG taste basis. For now, n,. will be a positive integer. We perform n + 1 RG 



blocking steps, labeled k = 0,1, ... ,n, following the blocking procedure of Ref. [15[. The 
special k = step is used to carry out the transition from the one-component to the taste 
basis, cf. Eq. (12. 4p . In this step the number of fermion degrees of freedom is not thinned out; 
in each subsequent step they are thinned out by a factor 2^ = 16. The partition function 
for this theory can be written as 

/n 
VUl[VV^'^ Br,{nr;U,{V^''^}) Det"" (Aa.te,n) • (3.1) 
k=0 

The notation here is as follows: The gauge fleld on the original lattice, with spacing a/, is 
denoted by U. The spacing of the k-th blocked lattice is = 2^^^af, and the gauge fleld on 
that lattice is V^''). The spacing of the flnal, coarse lattice is ac = 2""^^a/. The Boltzmann 
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weight for the collection of gauge fields, original and blocked, is -B„ (n,.; W, {V^^-*}). It 
is composed of three parts: the original gauge action, the gauge-field blocking kernels,^ 
and a short-distance contribution to the effective gauge-field action, n^.^S'e//, coming from 
integrating out the fermions on all lattices except the last one, where 



e 

k=0 



nDet(G,i). (3.2) 



The operators Dtaste,k and G^^ ^ are recursively given by 

DiLe,k = + Q^'^D;Xj.-iQ^'^^ > = 1, . . . , n , (3.3) 
G^' = Dtaste,k-i + «fcQ(')^g(') , A; = 1, . . . , n . 

The blocking parameter is of order l/ofc. The blocking kernel at the k-th. step, Q*^^^ = 
QCfe)^^")^^^-!))^ gauge-covariantly averages the fermion fields over 2^ hypercubes on the (A; — 1)- 
st lattice. For the A; = step, Dtaste,o = Dtaste is defined in Eq. f l2.4p . and Gq^ = Dgtag + 
olqQ^^^^ Q'^'^\ where = a and Q^^-' = Q are those introduced in Eq. (12.41) . Recall that the 
special = blocking kernel is unitary; all other blocking kernels are not. 

For small momenta, Q^'^^^Q^^^ ^ 1, and with ~ it thus follows that the eigenvalues 
of G^^ are at least of order l/flfc, making 5Seff a short- distance contribution to the effective 
gauge action.^ While this can be proved in the free case [34|, in the interacting case this is 
an assumption that is already necessary for the conventional RG picture to work in local, 



renormalizable theories. The nature of this assumption is discussed in detail in Ref. [15 
here we will assume it to be correct. It follows that 5Seff remains local when we take to 
be any real number.^ 

The fermionic contribution to long-distance physics then resides entirely in the n^-th 
power of the determinant of Dtaste,n in Eq. (13. ip . The problems with locality of the rooted 
theory originate with taking — 1/4 in this power. Our task will be to perform a faithful 
replica continuation at the level of the SET. As explained in the introduction, this is not 
straightforward. Calculations in the effective theories, the SET or ChPT, lead to explicit 
dependence on (for instance, through loops). But there is also implicit dependence 
through the couplings that appear in the effective theory, which is in general nonperturbative, 
and not known. 

Our strategy will be to first approximate the fourth-root theory by a local ( "re-weighted" ) 
theory. The fermions of this theory do not carry a taste degree of freedom; they are taste 
singlets. The multiplicity of taste-singlet fermions, n^, will always be chosen to match the 
fermion spectrum of the target continuum theory. Therefore we will never have to perform 
any "replica continuation" in n^; rather, Ug will always be kept a positive integer. In 
our construction, the unknown dependence of the couplings in the effective theory on the 
fermions will be due to the taste-singlet fermions only. 



We do not integrate over any of the gauge fields; this can be postponed to the end. The exphcit expression 



for B„ (n^;W, {V^^^^}) is given in Ref. [l5|. 



^ Much smaller eigenvalues are allowed, as long as the corresponding eigenmodes are localized on a distance 

of at most order a^. Such modes would not affect the long-distance physics. 
^ We keep rir in the range where the gauge coupling is asymptotically free. 
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The fourth-root theory will be reached from the taste-singlet theory by "turning on" 
the taste-breaking effects that introduce the nonlocal behavior. This is where a replica 
continuation away from the integers will be needed. Because of the smallness of the taste- 
breaking effects, the replica continuation will be under control. Indeed, we will show that 
to any order in oj, the dependence of the taste-breaking effects on is polynomial, with a 
degree less than the order in the a/-expansion. 

We start by splitting Dtaste,n into a taste-singlet part and a taste-breaking part with 
vanishing trace in taste space. 

An , (3.4) 
1 , 



■Dtaste,n 




Dinv,n 


Dinv.n 




1 








4 



D 



taste,n I 



where tr^^ denotes the trace in taste space, and 1 is the taste identity matrix. Following 



Ref. [15| we assume that, in the coarse-lattice theory, A„ scales like 



hAnW < ^ . (3.5) 

This estimate is valid modulo logarithmic corrections to the leading power-law scaling. For 
extensive discussions of this scaling assumption, we refer to Ref. [15|] (see also Refs. 



llip. Here we only observe that, in any theory with integer n^, this assumption is needed 
to establish that unrooted staggered fermions have the usually assumed continuum limit. 
However, by exploiting the proximity of the local re-weighted theory after a large number n 
of blocking steps, it was argued that the scaling (13. 5p is also valid in theories with fractional 
Uj.. In this paper, we will assume this to be the case. 

Using this split, we generalize the determinant in Eq. (13. ip to 

Det- (Aa..e,n) - Dct^ (A..,n) [n" ^ ' (3-6^) 

V / Det [Dinv,n> 



while also replacing 



Br. (n,;W,{V('-)}) ^ Bn (n,/4;W,{V('-)}) . (3.6b) 



The generalized theory reduces to Eq. (13. ip if we set = Ahj. and t = 1. This generalization 
has two important properties. First, if ng = and assumes physically interesting values, 
i.e., multiples of 1/4, then is an integer. Second, when n is large enough, D~^^A„ is 
small enough (in an ensemble-average sense) that we may expand 

""'Del^TAlf^^ = + ^A^nU^")] (3-7) 

-1) 



exp 



The parameter t interpolates between the taste-invariant operator Din^ n at t = and the 
(blocked) staggered operator at t = 1. In addition, t is a book-keeping device. The power 
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of t is, evidently, the same as the power of D^,^^^An. As we explain in detail in Sec. IIIIBI 
below, for the construction of the effective theories we may use the bound 

Pri,n^n||<^ = ^ = e„. (3.8) 

(The ~ sign has a meaning similar to that in Eq. (13. 5p .) We conclude that the t-expansion 
is an expansion in powers of a/ for the taste-breaking effects. 

For t = 0, the determinant ratio (13.71) collapses to one. The taste-invariant theory at 
t = is thus local for any integer Ug, and independent of n^. The staggered theory is 
reached by expanding as in Eq. (13. 7p . eventually setting t = 1. The rooted staggered theory 
is obtained by setting rir to a quarter-integer value. When we construct the SET to any 
finite order in a^, the maximal power of t will be limited by that order. ^ By Eq. (13. 7p . the 
maximal power of Ur is bounded by the power of t. (Because of taste-tracelessness of A„, 
the maximal power of is in fact strictly less than the power of t.) The maximal power of 
rir is thus (strictly) less than the order in a/. Therefore, at fixed Ug and to any finite order in 
ttf, the dependence of any correlation function on rir, and thus of the SET that reproduces 
it, will be polynomial. This implies that, at the level of the SET, the replica continuation 
in rij. to quarter-integer values will be well-defined, resulting in the "staggered SET with 
the replica rule." What this means is the following: We start with integer n^. The effective 
action is then given in terms of a set of Symanzik coefficients which are unknown functions 
of Us, but depend polynomially on rir (we may already set t = 1). With this action, one 
calculates correlation functions which again depend polynomially on Ur (to any finite order 
in a/), with rir dependence coming from the Symanzik coefficients and from loops. Finally, 
one sets rir = n-s/^, and the resulting correlation function is precisely that of the rooted 
staggered theory. The following subsections contain a more detailed argument on how this 
works. 

We comment in passing that, for t = 1, we may also interpolate between the taste-singlet 
local theory at rir = 0, and the (rooted) staggered theory at rir = ^s/i by varying rir 
instead of t. While the two ways of moving from the taste-singlet to the staggered theory 
are mathematically equivalent, we find the argument more transparent if the transition is 
done by varying t. 



A. The generalized theory 

In order to define the SET we first need a complete definition of the generalized staggered 
theory, coupled to sources in order to generate all correlation functions. Returning to integer 
rir, the theory defined by Eq. (13. 6p contains rig taste-singlet fermions with Dirac operator 
Dinv,n, nr generalized staggered fermions with Dirac operator Dinv,n + tAn, and irir ghosts 
with Dirac operator Dinv,n- Introducing sources H = {fj, rj, ff) and H = {r), fj, fj) for the taste- 
singlet, generalized staggered, and ghost fields respectively, we define the partition function 



^ Note that a/-dependence which does not involve taste-symmetry breaking may result from other sources 
besides the determinant ratio (|3.7[) . 
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of the generalized theory as 



/n 



(3.9) 



X Det"= D 



'mv,n 



Det""- ( A™,„ + tAr 



exp 



Det {Dinv,n) 



„ n 



fc=0 

X Det"= ( D,;„„.„ ) exp 



-rir 



X exp 
X exp 



Here / stands for the identity matrix, with dimensions as indicated by the subscript. This is 
a theory with two lattice parameters, Qc and af. Alternatively, we may trade af for the small 
parameter e„ of Eq. (13.81) . In the second expression we give the explicit expansions in the 
book-keeping parameter t. As explained above, for fixed Ug correlation functions expanded to 
some finite power in a/ are polynomial in n^. For t = 1, = ng/A, and 'fi = 'fj = fj = fj = 0, 
Eq. (13.91) is precisely the theory of degenerate, fourth-rooted staggered fermions. 

The generalized theory has a vector-like U{ns\4:nr) x U{nr) graded symmetry; U{ns\4:nr) 
acts on the taste-singlet and ghost fields, and U (rir) on the generalized staggered field. For 
t = the symmetry enlarges to f/(ns+4nr|4nr). The discrete symmetries include hypercubic 
rotations and axis reversal [3]. In the staggered sector, for t = 1 this is augmented by shift 
symmetry, and (softly broken) ?7(l)e symmetry in GWL form for each flavor. The vector 
and axial staggered symmetries expand to a U{nj.)e x U{nr)r chiral symmetry group (23| . 
There is no chiral symmetry in the taste-singlet and ghost sectors, because the GWL version 
of t/(l)e symmetry mixes the taste-invariant and noninvariant parts of the blocked staggered 
Dirac operator [l2|.^° 

We are now ready to discuss the SET for the generalized theory. As long as rir is a 
positive integer, the lattice theory is partially quenched but local, and we will assume that 
an SET for this theory exists in Euclidean space. The effective theory can be written in 
terms of continuum fields \1/ = {q,q,q) and \1/ = {q,q,q) for the taste-singlet, generalized 
staggered and ghost fields, respectively, as well as a continuum gluon field A^. As explained 
above, its parameters (the couplings multiplying each operator in the Symanzik expansion) 
are polynomials in Ur if we work to a finite order in aj; while their dependence on Ug is 



We remind the reader that Dinv,n and A„ in Eq. p.9p are defined in the RG taste basis, c/. Eq. (|2.4p . 
and not in the standard taste basis of Refs. (sl. Isst- 

It is sufficient to consider the SET in EucHdean space, since we will postpone the continuation to 
Minkowski space until after the continuum limit has been taken (11.] . 
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unknown. Only the Ug dependence survives in the continuum hmit, where the determinant 
ratio (13.71) collapses to one.^^ 

For general t, and Ur, the fundamental cutoff is the lattice spacing of the generalized 
theory, Oc- The SET is the effective theory for quarks and gluons with momenta much smaller 
than l/a^. However, the lattice theory contains an additional small parameter, e„ = af/ac, 
cf. Eq. fl3.8p . It will be useful for our purposes to think of the Symanzik expansion as 
an expansion in = e„ac, with Symanzik coefficients that depend on a^.^^ The effective 
theory can be divided into three different sectors, corresponding to three different types of 
operators that can occur. The (generalized) staggered sector consists of operators made out 
of staggered fields q and q only. Likewise, the taste-singlet-ghost sector consists of operators 

made out of the "auxiliary" fields \l/ = (g, q) and \1/ = (g, q) only. Finally there is the mixed 
sector, where each operator is made out of both staggered and auxiliary fields. (Of course, 
all operators may contain gluon fields.) 

In order to establish the validity of rSChPT in Sec. El we will not need to know the 
explicit form of the SET in full generality. In fact, we need only consider the staggered 
sector of the SET. Disregarding the auxiliary and mixed sectors, the resulting SET, defined 
in terms of the quark fields q and q and the gluon fields, is invariant under all symmetries 
of the generalized staggered operator Dinv,n + tAn- For t = this includes taste- replica 
symmetry U{4nr), while for t = 1 this includes the smaller group as well as softly broken 
U{l)e symmetry. 

For the remainder of this subsection we set t = 1, and thus Dinv,n + An = Dtaste,n reduces 
to the RG-blocked operator of Eq. (13. ip . Symmetries that act on the space-time coordinates 
often take a complicated form under RG blocking. In particular, shift symmetry is realized 
in a complicated way. First, the RG blocking leading to Eqs. (13. ip and (13.90 was started 
in the RG taste basis defined in Eq. (12.41) . and shift symmetry is thus realized as a gauge- 
covariant form of Eq. (I2.25p . Second, the transition to the RG taste basis was followed by 
n additional RG blocking steps. 

The physical consequences of any exact lattice symmetry of the underlying staggered 
theory, nevertheless, cannot be lost by RG blocking. The reason is the existence of a pull- 
back mapping of every coarse-lattice operator to a fine-lattice operator [l5|. For = Tis/4, 
where the taste-singlet and ghost determinants drop out, this mapping gives rise to exact 
equality of corresponding observables. In other words, the coarse-lattice observables are a 
subset of the original fine-lattice staggered observables. 

The pull-back mapping extends to 7^ Consider the expectation value of a 

product of coarse-lattice staggered fermion (and gauge) fields. By undoing the RG-blocking 
gaussian transformations of the fermions, this can be rewritten as an expectation value of 
a corresponding product of fine-lattice staggered fields (that depends in addition on the 
original and blocked gauge fields). Because the Boltzmann weight of the generalized theory 
contains the taste-singlet and ghost determinants, expectation values will not be the same 
as in the original staggered theory. But since the fine-lattice symmetries are unchanged, 
pulled-back coarse-lattice observables will still transform under all the staggered symmetries. 
Together with other observables constructed from the fine-lattice staggered fields, they must 
fall into representations of all these symmetries. This implies that the physical consequences 
of the full set of staggered symmetries remain intact. 



We observe that at nonzero Qc but a/ — > 0, i.e., in the hmit n — > 00, the lattice action is a perfect action. 
In the following subsection, we will argue that no negative powers of a/ can appear. 
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The t = 1 staggered-sector SET must therefore be invariant under all the symmetries 
listed in Sec. [Tll If we derive the SET using the taste basis some of these symmetries will 
take a complicated form. In particular, shift symmetry will mix different orders in a = af. 
But other continuum fields can always be chosen by suitable field redefinitions such that shift 
symmetry resumes the simple form of Eq. fl2.2p at the level of the SET. Moreover, a SET- 
level field redefinition will also eliminate any Oc-dependence of the SET that originates from 
the matching to the coarse-lattice interpolating fields. The only remaining dependence of 
the staggered-sector SET on Oc originates at this stage from the presence of the taste-singlet 
and ghost determinants in the underlying theory (13. 9p . 

Recall now that the group generated by the four elementary shifts 5^ contains translations 
by 2a f. At the level of the SET shift symmetry enlarges to the direct product of the group 
r4 and the continuous translation group. In the continuum limit af —>■ the discrete group 
r4 enlarges to the full taste/replica symmetry group SU^Arir) (with r4 embedded such that 
it acts identically on all replicas). 

The conclusion of the above arguments is that, for t = 1 and for any positive integer 
values of rig and n^,, the generalized staggered sector of the SET assumes exactly the same 
structure, as an expansion in the fine lattice spacing a/, as the standard staggered SET for 
Ur staggered fields. To order a'j, this SET is derived in Ref. jl^l (for Ur = 1) and Ref. jl^] (for 
arbitrary n,.), and is written down explicitly in Ref. jl^]. However, the Symanzik coefficients 
of the staggered-sector SET of the generalized theory are not the same functions of the 
parameters of the underlying theory as in the ordinary staggered SET. In the generalized 
theory, the Symanzik coefficients depend on and Oc, parameters not present in the ordinary 
staggered theory. Dependence on arises because of contributions from taste-singlet loops. 
In addition, the dependence (at fixed n^) of the Symanzik coefficients is different from 
that of the ordinary staggered SET, because of contributions from ghost loops. Indeed, the 
reason why the auxiliary sector was introduced in the first place, is that — unlike the original 
staggered theory — the SET of the generalized theory depends polynomially on rir to any 
order in a/, as long as is held fixed. 

We are now ready to make contact with the rooted theory. In order to reach the SET 
of the rooted theory we hold Ug fixed and choose t = 1. For any t, we may perform the 
replica continuation —>■ 71^/ A in any correlation function at any given order in the loop 
expansion. Indeed, because the Symanzik coefficients are polynomials in to any desired 
order in aj, this continuation from integer values of is well-defined. Now, recall that the 
taste-singlet and ghost sectors of the generalized theory (13. 9 p cancel (for vanishing sources) 
when we set Ur = Us/^- As explained above, this finally eliminates all the remaining depen- 
dence of the staggered-sector SET on the coarse spacing Oc, leaving only the dependence on 
the fine spacing a/. We have thus succeeded in constructing the replica-continued SET for 
the original blocked theory, Eq. (13. ip . for any quarter-integer value of n^, and to the desired 
order in a/. 

Putting everything together, we have shown that the familiar staggered SET for integer 



derived to order in Ref. j22|, |23|, and written down explicitly and extended to order 



Via the pull-back, the coarse-lattice operators may be regarded as a particular set of interpolating fields 
on the fine lattice as well. The freedom in making field redefinitions at the level of the SET thus parallels 
the freedom, discussed in Appendix B of Ref. jlO[, to choose different sets of interpolating fields on the 
fine lattice. 

^ For further discussion of the replica continuation, see Sec. IIII CI 
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a 



j in Ref. [2J], can be used to compute any correlation function of interest to the desired 
order in a/. The result should then be replica-continued to quarter-integer values of n^. 
This continuation provides the correct prescription for calculating any correlation function 
in the rooted theory from the staggered SET. Of course, in practice we will not know the 
precise coefficients of powers of in the Symanzik coefficients; indeed in practical situations 
the Symanzik coefficients must be treated as unknown numbers, to be fitted from numerical 
data. However, it suffices for our argument to know that the dependence is polynomial. 
When we continue in rir, we then need only continue the explicit Ur dependence coming 
from loops, giving a result as usual in terms of unknown Symanzik coefficients. 



B. Power counting 

A cornerstone in the argument of the previous section is the expansion in Eq. (13.91) . which 
is convergent if the norm of D^^^ ^ A„ is small enough. In this subsection, we consider this 
condition in more detail. There are two issues to be considered: the effect of insertions of 
An, as well as the size of the full object in which we expand, D^^^^An. 

In general, the SET for a lattice theory with lattice spacing a is constructed by matching 
correlation functions in an expansion in ap, with p <C 1/a a generic momentum, to the 
underlying lattice theory. To make the matching possible in perturbation theory, one should 
also take p ^ ^qcd- The Symanzik coefficients are extracted by computing suitable one- 
particle irreducible correlation functions in the lattice theory, taking all the (nonexceptional) 
external momenta to be of order p . For the part coming from the fermions, this amounts 
to expanding -D;~Jj around D'^^f, namely to an expansion in D'^nti^iatt — Dcont), where 
Dcont is the Dirac operator for the continuum- limit theory, and Diatt is the Dirac operator 
of the lattice theory. Because Diatt — Dcont is an irrelevant operator, we expect \\Diatt — 
DcontW ^ap^- Also, on dimensional grounds, ||-D^^J ~ l/p- Putting it together we conclude 
that \\D~^nt{Diatt — Dcont)\\ ~ o,p is the relevant estimate for the construction of the SET. 
Observe that this argument is insensitive to the long-distance physics, because the effective 
infrared cutoff on the loop momenta is p, and by assumption p ^ Aqcd- In particular, the 
estimates are independent of the quark masses. 

In the above argument we have implicitly assumed that the momentum flowing through 
a particular (sub-)diagram is of order p. This need not be true for sub-diagrams with a non- 
negative degree of divergence, where all ultraviolet momenta may contribute signiflcantly 
to the loop integrals. In general, counter terms will need to be added in order to absorb 
contributions from such diagrams; in a renormalizable theory there are only a flnite number 
of counter terms that need to be adjusted. Symmetries may exclude (some of) these counter 
terms. 

Let us now study how these general considerations enter the construction of the SET for 
the generalized theory (13.91) . Our starting point will be the t = taste-singlet theory. This 
theory is local, because Us is integer. In order to reach the generalized staggered theory from 
the taste-singlet theory, we have to expand the propagator {Dinv,n + tAn)~^ around D~^^^, 
and eventually set t = 1. The object in which we are expanding is thus Dj„^„A„. Since A„ 
is an irrelevant operator (c/. Eq. (13. 5p ). repeating the above general arguments leads to the 
estimate ||Dj^^„A„|| ~ Ufp, if the momentum flowing through the diagram is order p. 

As noted above, we must separately consider sub-diagrams with a non-negative degree 
of divergence. The contributions of such sub-diagrams depend crucially on the number of 
blocking steps n, as we now explain. 
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Consider first what happens for k = n = 0, namely, when we have performed only the 
first special RG step that takes the fermions from the one-component to the taste basis. We 
then have Oc = 2a f. When we extract the Symanzik coefficients from a lattice calculation, 
the loop momenta live on the coarse lattice. But since the coarse and fine lattice spacings 
differ only by a factor of two, the loop momentum can go as high as p ~ !/«/. In the 
divergent sub-diagrams we thus have „A„|| ~ 1. Indeed, for Oc = 2a/, the generalized 

staggered theory will develop 0(l/ac) = 0(l/a/) mass terms, since shift symmetry and 
U{l)e symmetry (for any t ^ 1) are broken at the (common) lattice scale. 

The situation is qualitatively different after a large number n of RG steps has been 
performed. Because the lattice calculation is performed on the coarse lattice, the maximal 
momentum that can flow through any sub-diagram is now of order l/oc, and one arrives at 
the estimate (13.51) for the magnitude of insertions of A„. The estimate \\D^^^ n^n\\ ^ ajp still 
holds, but, what has changed is that now the maximal value that p can reach is l/oc <^ l/c^/- 
The conclusion is that, for extracting the Symanzik coefficients, the appropriate estimate is 
just that of Eq. flXS]) : 

llAnlnAnll < a^/a, . (3.10) 

This estimate is valid in the taste-singlet, t = theory, on the same grounds as for any other 
local theory, and we will thus assume that it is valid nonperturbatively as well. This is all 
we need, because the staggered theory is constructed as an expansion in „A„ around 
the taste-singlet theory. 



We end this subsections with three comments. First, it should be noted that, in Ref. [15 
the bound 

llAinAnll < af/{ma',) (3.11) 

was used, with m the renormalized quark mass after n RG steps. Clearly, the bound (13. lip 
is far weaker than (I3.10p . and it implies that the chiral (m — > 0) limit can be taken only after 
the continuum (ay — > 0) limit. In Ref. [l^, this was necessary in order to place a uniform 
bound on the difference between any taste-singlet correlation function and the corresponding 
rooted correlation function on any (including the most infrared) scale, thereby establishing 
the existence of the (correct) continuum limit for the rooted theory. In contrast, assuming 
that the scaling (13. 5p holds, the bound (13. lip is much too generous for the derivation of 
the SET for the generalized theory (13.90 . as we have seen above. In particular, it follows 
that this SET is well-defined in the chiral limit, as is the chiral effective theory that can 
be derived from the SET. The requirement that the chiral limit for staggered fermions be 
taken after the continuum limit jl, 43, 44, 45, is then reproduced by calculations within 



staggered ChPT [45|. Note that, while Ref. [45j finds many standard quantities for which 



the limits commute in SChPT, other quantities for which the limits do not commute are 
also discussed. 

Our second comment is that the original staggered theory has no power divergences, 
because of shift and U{1)^ symmetry. This is therefore also true for the ra-times blocked 
staggered theory (13. ip . and for the corresponding SET. Moreover, for large n, the SET for 
the generalized theory (13.90 at arbitrary values for t G [0, 1) is related to the SET at t = 1 
by a convergent expansion in t, equivalently in e„ = af/ac- The implication is that, for all 



16 Thebrea; 
" See Ref. 



dng of shift symmetry is qualitatively the same as in the theory studied in Ref. [42 1. 



151 for a detailed discussion on how the coarse-lattice diagrammatic calculation is related to a 



calculation in the underlying fine-lattice staggered theory. 
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t, the SET for the generahzed theory (13. 9p has no power divergences in l/ct/, but only in 
l/oc- Examples of this are given in Sec. IIVI below. 

Finally, we remark that the framework introduced here resolves a concern, discussed in 
Ref. [lo|, about the renormalizability of the rooted staggered theory. The concern is the 
following: the complete notion of renormalizability requires not only that (infinite) countert- 
erms can be chosen to make amplitudes finite, but also that the finite parts of counterterms 
can be chosen to bring the theory into a given scheme. While we know that the staggered 
theory is renormalizable for integer tt.^, for non-integer this notion of renormalizability 
requires that the finite parts of counterterms, as well as the infinite parts, are polynomial in 
rir to any finite order in perturbation theory. In Ref. [13], the condition on the finite parts 
was introduced as an additional assumption, albeit a plausible one. Here, such a separate 
assumption is unnecessary. Under the assumptions of the RG approach [l5| , the taste-singlet 
(re- weighted) theory, defined by setting t = in Eq. (13.91) . is a local theory of Ug fermions, 
that moreover becomes a perfect-action lattice theory in the limit ^ 0, for any fixed ac- 
Thus one expects its renormalizability to follow straightforwardly by standard arguments. 
The rooted staggered theory is then reached by expanding in t, and setting t = 1 and 
Ur = Us/i. Because of the bound (I3.10p . the expansion in t just brings in positive powers 
of a/, and all finite (and infinite) parts of the counterterms are unaffected for any n,.. Thus 
the rooted staggered theory is renormalizable if the taste-singlet theory is. In addition, the 
two theories have the same counterterms. 



C. Partial quenching 

Unlike other lattice discretizations of QCD, the continuum limit of the rooted staggered 
theory is, inherently, a partially-quenched theory [lo|, [3, 0, This remains true when 



we consider the staggered sector of our generalized lattice theory (13.91) all by itself. Let us 
work out the example of a target theory with degenerate quarks. Our starting point is the 
generalized lattice theory with the same Ug, and with t = 1. In order to obtain the set of all 
correlation functions of the physical n^-flavor theory in the continuum limit, we need to let 
the combination of replica and taste indices of the external lines assume precisely Ug distinct 
values. This can, for example, be accomplished by fixing the taste index of the external 
legs to a single value (for example, 1), and letting the rephca indices take on Ug values 
(for example 1, 2, ... , Ug). Alternatively, we could use all four taste indices and only [ng/i] 
replica indices, where the square brackets denote rounding up to the next integer. (In this 
case, unless is already an integer, not all taste indices will be used in conjunction with 
each replica index.) Many other similar choices, as well as other types of embeddings for 
certain classes of physical correlation functions [l, [13], are also possible. Prior to the replica 
continuation, the lattice theory is local. The source term in Eq. (13. 9p must accommodate all 
the degrees of freedom, as specified above, that will be used in physical correlation functions. 
Therefore, we must consider only theories where n^, the (still integer!) number of staggered 
replicas, is not smaller than [ns/4]. 

When we perform the replica continuation we set the power of the staggered and ghost 
determinants in Eq. (13. 9p to rir = Ug/A. Since we have already set t = 1, if we turn off 
all sources, the partition function of the generalized theory reduces to the rooted partition 
function, in its RG-blocked dress (13. ip . During the replica continuation of any correlation 
function, by definition we hold fixed all indices of the external legs, including in particular 
the replica (and taste) indices. This means that the number of replicas in the source term of 
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Eq. (I3.9p must stay equal to or larger than [ns/4]. The mismatch created between the power 
of the staggered (or ghost) determinant and the multiplicity of the corresponding external 
sources means that the staggered sector has in itself been partially-quenched unless is a 
multiple of 4. 

After the replica continuation, the correlation functions of the EFT reproduce those of the 
rooted lattice theory to the same order in aj. We stress again that the replica continuation 
at the level of the EFT is well defined because, as we have shown, to any order in the 
rzr-dependence in the underlying lattice theory (I3.9p assumes the form of a finite-degree 
polynomial. 

In our above example, be it before or after the replica continuation, the replica x taste 
multiplicity of the staggered fields used to generated physical correlation functions is equal 
to or larger than 4[?t,s/4], which is to be compared with the physical flavors of the target 
theory. As a result, the total number of available valence degrees of freedom will in general 
exceed the physical number, and, when we finally take the continuum limit, the physical 
correlation functions will form a proper subset of the set of all (partially-quenched) correla- 
tion functions. This conclusion is in fact valid for any target theory. The only exception is 
a target theory in which the multiplicity of every mass-degenerate quark species is divisible 
by four, in which case the theory may be obtained in the continuum limit of an unrooted 
staggered theory. 

Another conclusion is that the partially-quenched representation obtained in the contin- 
uum limit is not unique. The only restriction is that the set of all partially-quenched corre- 
lation functions must be large enough to accommodate all the physical correlation functions 
of the target continuum theory. With the minimal choice of replicas on the external lines, 
[ns/4], the vector replica x taste symmetries are represented as a ?7(4[ns/4] | 4[n<j/4]— n^) 
graded group on the continuum- limit correlation functions. Had we initially allowed for 
n' > [^s/4] values of the replica index on the external legs, all the physical correlation 
functions of the target theory would still be reproduced once we performed the replica con- 
tinuation (followed by the continuum limit). But there would be more ways of embedding a 
given physical correlation function in the space of all correlation functions. Correspondingly, 
the replica x taste symmetries would be represented as an [/{An' \ An'—ris) graded group. The 
arbitrariness in picking a range n' > [ns/4] for the external-legs replica index thus entails 
the existence of infinitely many partially-quenched representations in the continuum limit, 
all of which share the same physical subspace. 

In the rooted theory, closed ("sea-quark") fermion loops as well as ( "valence- quark" ) 
fermion lines attached to external legs both originate from the same staggered fields. There- 
fore the sea and valence masses are equal, and there is no clear-cut distinction between 
the sea and valence sectors. This is a necessary condition for the emergence of a unitary, 
physical subspace in the continuum limit. 

In practice, it is often useful to explore unitarity-violating correlation functions in which 
the valence-quark mass is allowed to vary away from the sea-quark mass. This situation is 
what is usually referred to as partial quenching. As we have just explained, the continuum 
limit of the rooted theory is automatically a partially-quenched theory, albeit with equal sea 
and valence masses. If it is desired to study different sea and valence masses, it is straight- 
forward to add a (generalized-)staggered valence sector to the generating functional (13.91) . 



Correlation functions lying outside of the physical subset may exhibit various types of pathological 
behavior 16l. |46|. 
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by simply inserting a factor 

exp [rj, {Dl,^^ + t,A„)-i ® J„„) r],] (3.12) 

into the integrand. The superscript v on D^^^^ indicates that a different quark mass may 
have been chosen in the valence sector. For = 1 the valence sector has all staggered 
symmetries. Again, for a/ small enough, an expansion can be set up in t^, just as before. 

In Eq. (13.121) . rj^ and fj^ are sources for any desired number of valence (generalized) stag- 
gered fields. To avoid confusion we stress that, even if the valence-sector source term (I3.12p 
has been added to the generating functional (13.91) . we cannot dispose of the original source 
terms. The reason is that, if we want to consider the SET for both sea and valence quarks, 
we need sources for both in order to match the complete set of partially-quenched correla- 
tion functions between the lattice and the effective theory. With the valence sector (I3.12p 
in place, the replica x taste symmetries form an U{4n'+4ny \ An'+An^—ns) graded group in 
the continuum limit. (Of course, these symmetries will be softly broken by unequal sea and 
valence masses.) As before, n' is the number of distinct values of the replica index that we 
have allowed for the staggered fields with sea-quark mass on the external legs. 

In summary, we have seen that partial quenching occurs at three distinct levels. The 
generalized theory (13.91) is partially quenched to begin with, because, to keep the taste- 
breaking effects under control, we had to introduce a taste-singlet sector and a taste-invariant 
ghost sector. During the replica continuation, the staggered sector undergoes a second-stage 
partial quenching, created by the mismatch between the power of the determinant and the 
multiplicity of the sources. Last, if we are interested in different valence and sea masses, we 
need to introduce a "conventional" valence sector, cf. Eq. (I3.12p . 

IV. EXAMPLES 

It is instructive to consider some aspects of the SET to second order in aj in more detail. 
The SET can be written as an expansion in a/, t and n^, and thus takes the general form 

oo i J — 1 
1=0 j=0 k=0 

Here we already took into account that each power of t has to come with at least one power of 
a/, and that each power of has to be lower than the power of t (it cannot be equal because 
tTts{An) = 0). Equation (14.11) is manifestly polynomial in Uj. to any fixed, finite order in af. 
Here we allow all types of quarks (taste-singlet, generalized staggered and ghost) to appear 
on the external legs. The staggered sector is obtained by setting q = q = q = q = 0. The 
coefficients in Sij^k depend on both a^ and rig in all sectors. Because of this, one cannot in 
general conclude that terms linear in a/ have to be multiplied by dimension-five operators, 
etc. As already explained in Sec. IIII Al for t = 1 we may assume that a correlation function 
calculated in the SET does not depend on Oc if we set rir = ris/4: after the calculation. In 
this section, we will consider Ur integer. 



In this section we return to the theory defined by Eq. (|3.9p . The inclusion of valence quarks with a mass 
unequal to that of the sea quarks, as described in Sec. IIII C| is straightforward. 
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Because of the way the RG-blocked theory is constructed, for general t the preferred basis 
for (the generahzed staggered sector of) the SET is the RG-taste basis. Using this basis 
while restricting ourselves to the (generalized) staggered sector, and to i < 2, the expansion 
(14.1 p takes the explicit form 

S'^''^'{q,q,A;af,t,nr)= (4.2) 
+«/ [Sififiiq, q, A) + tS'i,i,o(g, q, A)] 

+a} [5*2,0,0(9,^,-4) +tS2,i,o{q,q,A) + 6*2,2,0 (g, g, ^) + nrt^5'2,2,i(g, g, ^)] • 

The nr-dependent term (the last term) is at this order the only one coming from the ex- 
pansion of the determinant ratio in Eq. (13.91) . The other t-dependent terms come from the 
expansion of the staggered source term in that equation. We note that S'j,o,o is taste invari- 
ant, because of taste invariance of the t = theory. Furthermore, >S'2,2,i is taste invariant too, 
because the factor of n^t^ originates from the determinant ratio in Eq. (13.91) . which does not 
affect the symmetry structure of the SET. The taste structure of the SET is determined by 
the external legs, which correspond to the source terms in Eq. (13. 9p . Since the two allowed 
insertions of A„ have been "used up" by the determinant ratio, only the taste-invariant part 
of the source term contributes to 5*2,2,1 • 

If we set t = 1 then, as discussed in Sec. IIII At there exist a field redefinition that brings 
gquad ^Yie familiar form of Ref. for rij. = 1, or to the form of Refs. [23|, |2J] for > 1. 



In particular, the redefinition removes the terms linear in a/. The Symanzik coefficients 



are equal to those of Refs. [23, l23|, if ^^^^ chooses Ug = Arir, a multiple of four 



For general rig and n^, the staggered SET is that of Refs. [23|, 12J], but the coefficients are 
different functions of n^.^" This form of 5''^""'^ is the one needed for the construction of 
rSChPT [2^, which we will discuss in Sec. El 

The taste-invariant operator Dinv,n has no chiral symmetry, even when the chiral limit 
is taken in the underlying staggered theory, and we would thus naively expect a linearly 
divergent mass term of the form gg/flc- However, for large n, the taste-invariant theory is 
close to the theory with t = 1 in the sense explained in Sec. IIII B[ In order to deviate from 
the t = 1 staggered theory, at least one power of a/, coming from an insertion of A„, is 
needed. Equivalently, the l/oc linear divergence has to be multiplied by at least one factor 
of e„ = af/ttc- In fact, even a mass term with magnitude ~ e^/oc = o//Oc cannot occur. To 
see this, note that we may write 

= D'^ — A + (A 3) 

taste,n 'inv,n inv,n ^ 'inv,n ' ' ' ' ' \ 

This shows that the order aj difference between the t = and t = 1 theories has to break 
taste, and therefore a taste-singlet difference has to be of order aj-. Singlet mass terms can 
thus only occur in 5*2,0,0 and 5*2,2,0, with opposite coefficients such that they cancel at t = 1. 

Next, let us consider nonsinglet mass terms, i.e., terms of the form gRq/a^ with some 
(momentum-independent) kernel K for which tTts{K) = 0. At order aj a nonsinglet mass 
term can only be part of 5i,i,o, because 5*1,0,0 is taste invariant. However, staggered sym- 
metries at t = 1 forbid such terms in 5*i,i,o, thus excluding this possibility. At order aj, a 



^'^ In particular, the Symanzik coefRcients of all taste-breaking four-fermion operators in the SET are 
independent of rir and depend only on n^. 
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nonsinglet mass term can only appear in tS'2,i,o + ^^82,2,0 because 5*2,0,0 and 5*2,2,1 are taste 
invariant. Let us assume that a bilinear qKq appears with coefficient ci in 5*2, 1,0, and with 
coefficient C2 in 5*2,2,0- Staggered symmetries then imply that tci + t^C2 = at t = 1, and 
thus Ci + C2 = 0. Any nonsinglet mass term at order is therefore proportional to t{t — 1). 
Simply put, there has to be a factor t in order to break taste symmetry, and a factor t — 1 
to break staggered symmetries, which include r4 and f/(l)e- 

In order to exclude various contributions to the nonsinglet mass terms in the above 
argument, we used the fact that mass terms cannot be introduced or removed by field 
redefinitions. As we now explain, the same is not true for operators of dimension five or 
higher: they cannot be excluded by arguments based on field redefinitions. With the taste 
basis of Eq. (14. 2p . we know from Eq. (12.281) that taste nonsinglet Wilson-like dimension-five 
operators will already appear in 5i,i,o. Of course, being nonsinglet, such terms will have 
to vanish at t = 0. In addition, because of staggered symmetries, a local field redefinition 
can be found removing such terms at t = 1. However, this same field redefinition applied 
to the SET at t 7^ 1 will, in general, introduce taste-breaking terms at t = 0. So, all we 
can conclude is that before the field redefinition such terms are proportional to t, while 
after the field redefinition they are proportional to t — 1 . We cannot conclude that they are 
proportional to t{t — 1). In the case of the mass terms discussed above, stronger conclusions 
are possible, because dimension-three terms cannot be removed by a field redefinition. 

V. STAGGERED CHIRAL PERTURBATION THEORY 

In this section, we will discuss the transition from the SET to staggered ChPT, or SChPT. 
For integer and = 4n.r the derivation was first given in Ref. j22| (for n,. = 1) and 
Ref. [2^ (for Hr > 1), and we refer to those papers for details on the explicit construction 
of the SChPT chiral Lagrangian. Here we will focus on the continuation to Ur = Ug/A, with 
rig as always a positive integer. 

A. The transition to staggered chiral perturbation theory 

In the previous section we explained how the appropriate SET for a rooted staggered 
theory can be constructed. Holding Ug fixed, the Symanzik coefficients are polynomials in 
Tir, and thus have no singularities at quarter- integer values of n^. The rooted staggered SET 
is obtained as a replica rule: calculate correlation functions to a given order in aj, then 
set Hr = ns/4:. For the next step — the transition to ChPT — we must again retain both 
and fir as independent variables. In ChPT, as for the SET, the replica continuation in rir 
will be well-defined at fixed n^, and SChPT with the replica rule, namely rSChPT, will be 
recovered after the continuation to Ur = rig /A. 

When we calculate correlation functions using the SET for the generalized theory (13.91) . 
dependence on occurs in two ways: through the polynomial dependence of the Symanzik 
coefficients, and through fermion loops. Once we have calculated a certain correlation func- 
tion to some order in aj and to a given order in the loop expansion, the dependence on n^. is 
thus explicitly known. Technically, this dependence will not be a polynomial, because only 
the inverse quark propagators, and not the quark propagators themselves, depend polyno- 
mially on Ur- However, each quark propagator can be re-expanded around that of the t = 
theory in terms of a/, and thus n,-, just as in the underlying lattice theory (Eq. (13. 9p ). 
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This sets the stage for the derivation of the appropriate chiral theory for QCD with rooted 
staggered fermions. The continuum chiral theory is an effective theory for low-energy scales 
where only Goldstone bosons can appear on the external lines. It can be organized as an 
expansion in p/Ay, where ~ 1 GeV is the chiral scale separating other hadrons from the 
Goldstone bosons |47| . The chiral effective theory can be generalized to include discretization 
errors, in an expansion in a = a/. The chiral effective theory is to be constructed by matching 
its correlation functions to those of the underlying theory in a double expansion in p/A^ and 
ajp. In practice, the low-energy constants (LECs) of the chiral theory cannot be calculated 
by analytic methods, and are determined by fitting experimental or numerical data. 

For positive integer Us and n^, the underlying lattice theory is local, as is the SET, and 
the transition to the chiral theory is more or less standard [22, [23, 48, In addition, the 

estimate fl3.10p is still expected to hold, even though it cannot be checked in perturbation 
theory, because in this case the correct degrees of freedom for p ^ A^ are no longer quarks 
and gluons. Using the expansion (13.91) just as in Sec. IIIH this implies that the LECs of the 
chiral theory again have to be polynomials in n^. Finally, setting t = 1 and performing 
the continuation to = 'n.s/4 we recover the replica-continued SChPT, or rSChPT, of 
Refs. [i3,[23. 

We assume here that the contributions of ghosts and taste-singlet quarks in the sea will 
cancel to all orders in the partially quenched ChPT once we put = ns/4. All differences 



between the current rSChPT and the standard rSChPT [ij, |23|| (which does not have the 
taste-singlet and ghost sectors) will then disappear in the limit = ns/4, as long as we 
choose not to put ghosts and taste-singlet quarks on the external lines. Since the ghost 
and taste-singlet Dirac operators and masses are identical, this cancellation is trivial at the 
QCD level, but not completely trivial beyond one loop at the chiral level. We believe, 
though, that the cancellation is almost certainly true order by order in SChPT, and that 
it will probably be possible to construct a "quark flow" proof of this. This completes our 
argument that rSChPT is the correct chiral theory for QCD with rooted staggered fermions. 



B. An example 

It is instructive to see how our approach works in a concrete example. We will re-consider 
the leading-order contribution in rSChPT to the connected scalar two-point function, previ- 
ously described in detail in Sec. 6 of Ref. [13] ■ Adding a scalar source s{x) to the generating 
functional, this two-point function is defined as the connected part of the second derivative 
with respect to this source (setting s = after taking the derivatives). Adapting it to our 



generalized theory, Eq. (27) of Ref. [ij] takes the form' 



,23 



/ riLi ^^V(^) {^) Det"^ ( Aa.te,n + s®l) Det^"^-^"") (d„ + 

Z{s) = j-^ ^ , (5.1) 

JVUULil^VO'^ Bn (t) Det"" (Aaste,n)Det("-^"'-) (A™ 



Again, the only element of this transition that is not absolutely standard is the assumption that all steps 
can be carried out for partially-quenched theories, since the generalized theory ()3.9|) is partially quenched. 
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We thank S. Sharpe for emphasizing this point to us. 



^•^ The connection with the method and notation of Ref. 14| is explained in Sec. IV Cl 
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where we only indicated the dependence of Bn exphcitly, cf. Eq. (13. 9p . Here we have 
chosen t = 1, but have not yet set Ur = ng/^. It is important to keep Ur integral at this 
stage in order to develop the chiral theory; keeping n.^ 7^ Ug/^ also allows us to highlight the 
different ways in which rig and appear. 

In Eq. f l5.ll) . we are starting from the fact that correlation functions generated in the 
rooted staggered theory by the taste-singlet meson source s{x) ® 1 are identical, in the 
continuum limit, to the desired correlations generated by s{x) in the target QCD theory. 
(See Eq. (12) of Ref. (sj.) Note, however, that we have coupled s{x) not only to the staggered 
quarks but also to the ghost and taste-singlet quarks. This keeps the expansion in under 
control because the staggered and ghost contributions differ only by the small taste- violating 
term A„. Requiring that the taste-singlet and ghost quarks cancel at rtr = Ug/A then implies 
that s{x) also couples to the taste-singlet quarks. 

Even without a replica continuation, the lattice theory defined by Eq. (15.11) is, as we 
discussed already above, a partially quenched theory with rir staggered fermions, Ug taste- 
singlet fermions, and An^ taste-singlet ghosts. It differs from Eq. (13.91) in the way it is coupled 
to sources. Of course, the correlation functions that are generated by taking derivatives with 
respect to s{x) can also be generated by taking joint derivatives with respect to H{x) and 
H{x) (with one each for each space-time point). Regardless of which type of source is 
used, the dynamics is that of the sea-quark loops, and is controlled by the determinants in 
Eq. (13. 9p . Since in this subsection we are only interested in the scalar two-point function, the 
formulation with the source s{x) is simpler. Note that here we need the complete effective 
theory, including taste-singlet and mixed sectors, because the source s{x) couples to all 
quarks. 

At leading order in ChPT, the scalar two-point function consists of a sum over one- 
loop diagrams, with pseudo-scalar mesons on the loop {cf. Fig. 2 of Ref. [3l)- Since s{x) 
couples to all bilinears, staggered, taste-singlet, and ghost, all types of pseudo-scalar mesons 
contribute to these diagrams, including fermionic mesons made out of quarks and ghosts, 
and mesons made only out of ghosts. Because the taste-singlet quarks and ghost have the 
same Dirac operator Dinv,n, the result for the scalar two-point function that we will give 
below is that of a theory with Ug — An^. taste-singlet quarks, irrespective of the value (and 
in particular, sign) of ng — Arir- In the interest of brevity, therefore, the discussion below 
will simply assume that we are dealing with a theory with a positive number ng — An^ of 
taste-singlet quarks (as well as staggered quarks). 

In Ref. [13] it was shown that, as expected, in the one-flavor theory (for which Ug = Aur = 
1) only the non-Goldstone, heavy pseudo-scalar taste-singlet state (the "77"') contributes to 
this two-point function in the continuum limit, despite the presence of flfteen additional 
light pions in the underlying staggered theory. That this has to happen follows from the 
general discussion given in Ref. [S]- Here we will not repeat the details of the calculation 



given in Ref. [IJ] , but only keep track of how the results change in the generalized setup of 



the present paper, and see how and Ug appear in the flnal result. With Ref. [ij], we keep 
the singlet pseudo-scalar state in the calculation for pedagogical reasons. 

There are now three kinds of pions, those made out of staggered quarks, those made out 
of taste-singlet quarks, and "mixed pions," made out of staggered and taste-singlet quarks. 
The leading-order masses of the pseudo-scalars in the staggered sector are given by 

Ml = 2/im + ajAs , (5.2) 

where S G {I,C,fj.,iC,fiC,u{^J' > ^),i(,fj.^5,(,5} labels the taste of each of the sixteen staggered 
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pseudo-scalars (for each replica), and the Ah are four LECs^'' representing the taste spht- 
tings; m is the quark mass. Then there are pions made out of only taste-singlet quarks, with 



9^ 

mass 



M/, = 2/im + ajAts . (5.3) 
Finally, there are mixed pseudo-scalars made out of one taste-singlet and one staggered 



quark. The mass of the latter can be parametrized, to leading order, as (50 



M^.^ = 2/im + a^fArnix , (5.4) 



with, in general, Amix 7^ Ats- The fact that the mass of the mixed mesons does not depend 
on their staggered taste follows, as in Ref. HO], from shift symmetry, which forbids taste- 



violating staggered bilinears, and therefore forbids taste- violating four-quark operators with 
one staggered and one taste-singlet bilinear. Note that all the above masses (in particular. 
Mi) are the pseudo-scalar masses before including the effect of the anomaly. 

The LECs fi, Ah, A^^ and Amix have unknown dependence on Ug, but do not depend on 
rir. For /i this is obvious, because it is a continuum LEG, and the continuum theory does not 
depend on n,. at all, but only on n^. (Recall that, in the continuum limit, the determinants 
ratio (13. 7p goes to one.) Because Ah represents an order aj effect, it can, according to our 
general arguments, be at most linear in rir. In practice, it is independent of n^, because 
symmetry-breaking terms of order a'j in the SET do not originate from the determinant 
ratio but only from the source term in Eq. (13.91) (c/. the discussion below Eq. (14.21) ): similar 
arguments apply for Amix and A^^. At higher order there will be r;,j.-dependent corrections 
to Eqs. (15. 2p through (15. 4p coming from insertions of the operator 6*2.2,1 in Eq- (14.21) . The 
taste-singlet and mixed mesons also contribute to our scalar two-point function as long as 
rir 7^ n^/4. 

Of course, the singlet pseudo-scalar (the "r/"') will not be a Goldstone boson. It will pick 
up a mass that does not vanish in the chiral and continuum limits. In the continuum limit, 
the rj' mass is given by 

M2=2/xm + n,^, (5.5) 

where ml is the double-hairpin parameter (c/. Ref. (23|). 

Again, since the continuum limit does not depend on rir, the parameter rriQ does not 
depend on r?,,,.^^ Away from the continuum limit, mixing takes place in the neutral meson 
sector because of different scaling violations in Mf and M^^. This mixing leads to the 



Ajj; = because this taste corresponds to the exact Goldstone bosons. 

The operator Dinv,n has no chiral symmetry, and the taste-singlet quark mass is additively renormalized 
by an amount of order aj (see Sec. lIVp . The quantity Ats represents the effect of this renormalization on 
the meson mass. In the case of the mixed pseudo-scalar mass, Eq. (|5.4p . such renormalization is absorbed 
in Amix, which must be present in any case. 

There are in general corrections of order aj, as well as momentum-dependent contributions, to this 
parameter, but they do not invalidate our conclusions. Following Ref. [14ii] . other hairpin contributions of 
order aj will be ignored as well. 
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appearance of pseudo-scalar mesons with masses M± given by 



M: 



II ml 



2\ 2 



m, 



_2 



(5.6) 



In the continuum hmit, A = and M| = Ml = 2fim, so the expression for simphfies 
to Eq. fl53|) . 

In order to give the expression for the scalar two-point function, we define single-particle 
propagators 

1 

A = ts, mix , 



DAip) 

and hairpin "double poles" 



p2 + M2 ' 



XiAp) 

Xts,tsip) 
Xl,ts{p) = Xts,l{p) 



1 



p2 + M^ 



{p"^ + Ml){p'^ + Ml) p"^ + Mf ' 

1 + M| 

(p2 + M2)(p2 + p2 + ' 

1 

(p2 + M2)(p2 + M2) ■ 



(5.7) 



(5.8) 



For A = 0, all hairpin double poles become equal, and Dj{p) = Dts{p)- 

The result for the Fourier transform G{p) of the scalar two-point function is 



G{q) = /i^ ' 



(5.9) 



+16nr{ns - Anr)Dmix{p)Dmix{p + q) + 2(ra^ - AnrfDts{p)Dts{p + g) 



-8nr^ {Djip)Xjj{p + q) + Dj{p + q)Xjj{p)) 



-2{ns - Aur)^ {Dtsip)Xts,tsip + q) + Dts{p + q)Xu,ts{p)) 



rrir 



2\ 2 



32nlXjj{p)Xjj{p + q) + 2(n, - 4n,)2Xi,,t,(p)Xi,,t,(p + q) 



+ IGririris - Anr)Xi^ts{p)Xi^ts{p + g) 
The explicit factors rn^/?) can be eliminated from this expression by using the relation 

.2 1 



m: 



^ = _ (M^ + m! - - Ml) . 
3 n. 



(5.10) 



As discussed above, if we expand out the masses Ml in powers of a^, the rir dependence of 
Eq. (15. 9p is polynomial. The dependence is not polynomial because the LECs /x. As, A^^ 
and Amii depend on Us implicitly in an unknown way. 
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Let us compare the result (15.91) to a similar calculation, done in the taste-singlet theory 
obtained by replacing Dtaste,n with Dinv,n in Eq. (15.11) . To order aj, this corresponds to 
setting M| = M^^^ = M^^. The expression for (c/. (5.6)) again simplifies to (5.5), 
except that 2/im is replaced with Mf'g, because M^^ may still include discretization errors. 
Instead of Eq. (15.91) we now arrive at 

1 1 

Ml,s = Ml + n,^. 

As expected, this result is n^-independent. The first term on the right-hand side is recognized 
as the anticipated contribution of the — 1 degenerate Goldstone pions of a theory with 
Us (mass-degenerate) fiavors. 

Replacing Dtaste,n with Dinv,n means that the product of determinants in the denominator 
of Eq. (15. ip collapses to Det"''(Dj„t,_„), with a similar simplification in the numerator. Our 
calculation thus explicitly demonstrates how we may consider the rooted staggered theory 
as a local taste-singlet theory with small, nonlocal corrections of order a^, which, to any 
fixed order in a/, are polynomial in n^.^^ Our example also illustrates how the nonlocality of 
the rooted staggered theory manifests itself in the low-energy EFT: while Eq. (15.1ip satisfies 
unitarity, Eq. (15. 9p . at 7^ 0, does not. This is most easily seen by noting the presence 
of the minus signs multiplying various terms in Eq. (15.90 . in what should be (in a unitary 
theory) a positive definite correlation function. 




C. Comparison with Reference [14] 

The present work may be compared with the complementary argument for the validity 



of rSChPT given in Ref. [IJ]. That argument starts from ChPT for a rooted theory with 
four degenerate fiavors of staggered fermions, which thus describes four mass-degenerate 
quark species. The underlying lattice theory is local, trivially, because it contains the fourth 
power of the fourth- rooted staggered determinant. Staying entirely within the ChPT frame- 



work, Ref. [1J| then treats the nondegenerate case by perturbing in the quark masses. An 
assumption of the analyticity of the expansion around positive quark mass is required at 
this point. In addition, the replica rule (called the "replica trick" in Ref. jl3]) needs to be 
introduced because the theory becomes nonlocal as one moves away from the degenerate 
limit. Finally, one of the four masses can be made so large that that quark decouples from 
the chiral effective theory (at which point it can be thought of as the charm quark). Using 
an assumption about the details of decoupling, one arrives at rSChPT for three light quarks. 
The decoupling assumption leaves a small potential loophole in the argument of Ref. [l^ 



By making use of the general sources in Eq. p.9p this conclusion applies to any physical correlation 
function of interest. A by-product is that the generalized theory (|3.9p provides a alternative framework 
to that discussed in Appendix B of Ref. [l^] for solving the "valence rooting" problem. . 
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While the three-flavor chiral theory goes over, in the continuum hmit, to the standard three- 
flavor chiral theory of QCD, it is not guaranteed that the LECs have the same numerical 
values as in QCD. (In the initial four-flavor case, the correctness of the LECs is guaranteed, 
however.) 

Here, we have started instead from the fundamental lattice theory (in RG-blocked form) 
and have shown how rSChPT may be derived from it, via the SET. The replica rule is given 
deflnite meaning in the fundamental theory, so its appearance in the EFTs is completely 



natural. In contrast, the replica rule in Ref. [IJ] has, by construction, meaning only at 
the chiral level. It is for that reason that a distinction was made in Ref. ij] between the 
power of the staggered determinant at the QCD level, which was called R, and the number 
of replicas introduced at the chiral level, rij.. Here, because the replica rule is justifled at 
the QCD level, we need make no such distinction. We do however need to introduce the 
number of flavors of the taste-singlet quarks, Ug, which affects LECs in a nonperturbative 
(and hence unknown) way, in order that the dependence be completely controlled (indeed, 
polynomial). Thus Ref. [l^ and the current work represent two different generalizations of 
the staggered theory. In the limit R = rir = Ug/A, the two generalizations agree. Since this 
is the limit we need to take at the end of any rSChPT calculation, it is clear that the two 
versions of rSChPT give the same results. 

Another advantage of the present approach is that it allows us to dispense with the 
assumptions about decoupling and about the analyticity of the mass expansion. This means 
that the current argument closes the loophole mentioned above. The continuum low-energy 
constants are automatically those of QCD with the correct number of flavors. 



On the other hand, the current argument, based as it is on Ref. [15|, inherits the assump- 
tions of that work. The key assumptions have already been mentioned in the Introduction 
and explained in Sec. IIIII They are that: 

• The effective action SSejf, generated by integrating out fermions on flner lattices, is 
local. 

• The perturbative scaling laws apply, implying that the dimension-flve taste-breaking 
operator A„ goes to zero like aj (times logarithms) in the continuum limit. This in 
turn is based on the highly plausible assumption that the theory is renormalizable to 
all orders in perturbation theory for any n^. 



The assumption of taste-symmetry restoration is needed in Ref. [14| too, but only for 
integer n^, where the scaling argument is completely standard. The argument of Ref. [3] 
works entirely within the chiral theory, and the resulting rSChPT then implies the symmetry 
restoration (in the chiral sector) for the rooted case. We also note that, in the RG framework, 
there is an alternative route to establish the validity of the continuum limit while relying 
only on the scaling of A„ in the taste-singlet (re- weighted) theory Since the latter 

theory is local by the flrst assumption, the validity of the scaling assumption needed for the 
RG treatment is very plausible. We remind the reader that there is considerable numerical 



evidence for the continuum restoration of taste symmetry in the rooted case [13|, ll8|, ll9|, |2 



21 . 



We again are assuming that the contributions of ghosts and taste-singlet quarks in the sea cancel to all 
orders in partially quenched ChPT once there are the same number of ghosts and taste-singlet quarks, 
i.e., once = ns/4. 
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Both the present arguments and those of Ref. rely heavily on the validity of the 



standard partially quenched chiral theory [16|] for describing partially quenched fundamental 
theories that are local. We also need to assume here that the SET exists for partially- 
quenched theories, as long as the lattice theory is local. 

The calculation of the scalar two-point function, presented in Sec. IV Bt may now be 



compared to the corresponding calculation in Sec. 6 of Ref. [Ij]. Note that Ref. considers 
only the one-flavor case as an example, so to make the connection, we must put = 1. 
The result here, Eq. (15.91) . then corresponds directly to Eq. (41) of Ref. [3]. We can in fact 
make the connection at the quark flow level: The first two lines of Eq. (15.91) correspond to 
Figs. 3(a) and (d) of Ref. [ij], the next two lines correspond to Figs. 3(b) and (c), and 
the last two lines correspond to Fig. 3(e). It is straightforward to check that, if we set 
Ur = Us/ A = 1/4: in Eq. fl5^ . and = = 1/4 in Eq. (41) of Ref. [li], the results are 
identical. 



VI. CONCLUSION 

In this paper we presented a theoretical argument that rSChPT 2^ is the correct chiral 
theory for QCD with rooted staggered fermions. Much evidence in favor of this claim alr eady 



existed, both on the theoretical side [Ij] , as well as on the numerical side [13|, ll8|, Il9|, 120 
Here we showed that it is possible to extend the usual construction of the Symanzik effective 
theory and chiral perturbation theory, to the rooted staggered case. Our arguments apply 
equally well to any staggered quark action that has the usual staggered symmetries, for 
example standard (unimproved) staggered pj], Asqtad ^], HYP [29|], Fat7bar |30j, or HISQ 
[m quarks. The version of staggered quarks used will not effect the form of the discretization 
effects summarized by the effective theory, but does effect the size of these effects, which is 
reflected in the size of the LECs. 

The effective theories are first constructed for a taste-singlet local theory with physical 
fermion flavors (the t = theory of Eq. (13.91) ). The rooted, nonlocal staggered theory is 
then reconstructed as an expansion in the lattice spacing of the underlying staggered theory 
{i.e., aj), by moving smoothly from t = to t = 1. In this framework, the dependence on 
is polynomial to any finite order in a/ and to any finite order in the loop expansion. The 
effective theories, however, are in the first instance only known at integer values oirir, where 
they are fairly standard. The polynomial dependence on rir allows us to to make the replica 
continuation of any correlation function, computed order-by-order in the effective theory 
for integer n^, to rij. = rig/ 4. Once the value = Ug/A is reached, the correct correlation 
functions of the underlying rooted lattice theory are recovered. 

The ability to extend standard techniques for the derivation of the SET and ChPT 
to rooted staggered fermions does not preclude various sicknesses in the rooted theory at 
nonzero a/. Indeed, in Ref. 12] we argued that the rooted theory is nonlocal at nonzero a/, 
due to the taste-breaking induced splittings in hadron taste multiplets. It is essential that the 
replica-continued SET and SChPT reproduce the nonlocal behavior. This happens because 
loop corrections calculated in these theories have to be continued to noninteger number of 
staggered replicas as well, and the replica-continued amplitudes cannot be reproduced from 



For the SET, the relevant loop expansion is the one in fermion loops; for ChPT it is the chiral loop 
expansion. 
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any local Lagrangian. An explicit example of this was worked out in Sec. 6 of Ref. [IJ]; we 
revisited this example in Sec. IVBl in our generalized framework. 

It is important to list the assumptions that underlie our arguments. The most important 
assumption is that QCD with rooted staggered fermions has the desired continuum limit. 
This conclusion, in turn, is based on a number of technical and testable assumptions, as 



explained in detail in Ref. [15| (see also Refs. [10|, lUi)- If this conclusion were to turn out 
to be incorrect, that would also invalidate the analysis presented here. Turning this around, 
we consider the success of fitting high-precision numerical results with rSChPT as direct 



evidence that the conclusion of Ref. [15[ is, in fact, valid. 

In order to keep the replica continuation under control, in Eq. (13. 9p we temporarily 
treated the number of dynamical quarks in the theory (n^) and the power of the staggered 
determinant (n^) as independent. Because ghosts are needed, we also have to assume 
that the construction of the SET and ChPT goes through in the standard way for partially- 
quenched (but local) theories. This second assumption is very common in applications of 
EFTs to lattice QCD. However, one should keep in mind that, while partially quenched 
ChPT [l6[ is by now standard, its foundations are not as firm as for ordinary, unquenched, 
ChPT. See Ref. j5l| for a discussion of this point. 

A third assumption is the technical observation that D~^^ „A„ has to scale as a/p, with p 
the momentum scale at which a correlation function in the effective theory is matched to the 
underlying theory. An exception are short- distance contributions coming from sub-diagrams 
with non-negative degree of divergence in which D~^^A„ can become as large as aj/uc at 
most. The end result, the estimate (I3.10p . is crucial for establishing that the n^-dependence 
of the generalized theory (13. 9p is polynomial, to any finite order in a/.^° Again, we consider 
this assumption as noncontroversial, because it underlies the standard derivation of EFTs 
for local lattice theories, and because it is used only in the t = theory, which is local by our 
first assumption. The weaker, quark-mass dependent bound on D~^^.^An used in Ref. [ll] 
is not needed for the derivation of the effective theories, and both the SET and the chiral 
theory are valid in the chiral limit. We emphasize here that the physically sensible approach 
for any staggered theory (rooted or not) is to avoid the region m <C OjAq^^, where lattice 



artifacts may dominate [8|, |44, |45 



In the actual construction of a SET or a chiral theory, use is made of the symmetries of 
the underlying theory. Particularly important symmetries for staggered fermions are U{1)^ 
chiral symmetry and shift symmetry, and we discussed in detail how these are realized at 
the level of the SET. Generalizing a result previously derived to order a'j in Ref. [l^l, we 
showed that for the SET, shift symmetry enlarges to the direct product of the continuum 
translation group and the finite discrete group Since this observation holds for the SET, 
it also holds for any EFT derived from the SET. Finally, we note that our arguments also 
apply to the cases of rSChPT with baryons or heavy-light mesons. 
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